Prom multiferroics to cosmology: Scaling behaviour and beyond in the hexagonal 

manganites 
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We show that the improper ferroelectric phase transition in the multiferroic hexagonal manganites 
displays the same symmetry-brealcing characteristics as those proposed in early-universe theories. 
We present an analysis of the Kibble-Zurek theory of topological defect formation applied to the 
hexagonal manganites, discuss the conditions determining the range of cooling rates in which Kibble- 
Zurek behavior is expected, and show that recent literature data^ are consistent with our predictions. 
We explore experimentally for the first time to our knowledge the cross-over out of the Kibble-Zurek 
regime^ and find a surprising "anti-Kibble-Zurek" behavior. 



A key open question in cosmology is whether the 
vacuum contains topological defects such as cosmic 
stringsS"— . The formation of topological defects during 
phase transitions in the early universe - believed in grand 
unified theories to occur within around 10~^* s after 
the Big Bang - was proposed by Kibble^, who derived 
the symmetry requirements for formation of topological 
defects in a physical manifold. In systems where such 
defects are symmetry allowed, their density can be es- 
timated using the Zurek scenario^, which uses causal- 
ity arguments to develop scaling laws for the density 
of defects formed as a function of the rate of quench- 
ing across the phase transition. The resulting combina- 
tion of symmetry requirements and scaling laws is termed 
the Kibble-Zurek mechanism, and in principle should de- 
scribe a phase transition in any system with the required 
symmetry properties, provided that other effects do not 
dominate the kinetics of topological defect formation. 

Attempts to demonstrate Kibble-Zurek scaling in 
condensed-matter systems to date have proved chal- 
lenging, however, and the "ideal Kibble-Zurek system" 
has remained elusive. Zurek's original paper- discussed 
the analogue between cosmic strings and the vortex 
cores formed in a quench-induced phase transition from 
normal-state to superfluid '*He. However the correspond- 
ing experimentiSt^^ yielded large deviations from the pre- 
dicted behavior, probably because of thermal effectsi^. 
In ^He the symmetry breaking is closer to that pos- 
tulated for the early universei^ii^, but the density of 
topological defects can only be inferred indirectly and 
many assumptions must be made to compare with pre- 
dictions. In superconducting Nb rings the density of vor- 
tex cores in the superconducting current led to a differ- 
ent scaling exponent than that predicted by the Kibble- 
Zurek mechanism - again, experimental artifacts were 
held responsibleiS, - and no flux-line formation at all 
was observed in experiments on high-Tc superconducting 
thin filmai^. Bose- Einstein condensates could in principle 
provide a suitable system but are so far subject to experi- 
mental limitations^^. Perhaps the best candidates to date 
have been liquid crystals, in which the diffraction pattern 
formation in nonlinear-optic experiments has been shown 



to exhibit power-law scalingi^ and promising studies of 
defect dynamics have been performed^^, but strong inter- 
actions between the defects continue to cause difficulties. 

Here we propose the multiferroic hexagonal mag- 
nanites, i?Mn03, R = Sc, Y, Dy - Lu, as a model 
system for testing the Kibble-Zurek mechanism. The 
hexagonal manganites have attracted interest because 
of their unusual geometrically-driven improper ferroelec- 
tricity, which allows for the simultaneous occurrence of 
magnetic ordering22i2i, as well as unusual couplings^^"— 
and functionalities^^ at their domain walls. In this work 
we show that the unusual nature of the improper geomet- 
ric ferroelectric phase transition also sets both the correct 
symmetry conditions for Kibble-Zurek behavior to man- 
ifest, as well as the physical properties for readily detect- 
ing it. In addition, the relevant time-, temperature- and 
length-scales fall into a range that allows exploration of 
the Kibble-Zurek regime, as well as the crossover out of 
it. 



I. SYMMETRY AND PHYSICAL PROPERTIES 

OF i?MN03 

First we describe the properties of i?Mn03 that are 
relevant for testing the Kibble-Zurek mechanism, par- 
ticularly the symmetry properties of the phase transi- 
tion. The structure of i^MnOa consists of planes of 
MnOs trigonal bipyramids separated by planes of R ions 
which form a hexagonal mesh (Fig. l(a))^. In the high- 
temperature paraelectric phase, the space group is cen- 
trosymmetric PG^/mmc. At Tc ^ 1400 K (the exact 
value depends on the R ion) a spontaneous symmetry 
breaking occurs, with the condensation of primarily two 
phonon modes with distinct irreducible representations 
of the high-symmetry structureSi^. First, a mode of 
symmetry, which involves a trimerizing tilt of the trigo- 
nal bipyramids and is the primary order parameter (Fig. 
1(b)). Since the if 3 mode can condense about three dif- 
ferent origins, and the tilt can be in the "in" or "out" 
direction, six trimerization domains are formed; these 
have been shown using high resolution transmission elec- 





FIG. 1: (a) High-symmetry PG^/mmc structure of RMnOs 
before the onset of trimerization. (b) Action of the K3 trimer- 
ization mode on the R ions and MnOs trigonal bipyramids. 
The insets below the main structure are to emphasize that 
outward trimerization results in a downward shift of the cor- 
responding R ion, whereas inward trimerization results in an 
upward shift, (c) The subsequent additional displacements of 
the R ions (blue arrows) in the mode provide the ferro- 
electricity. Note that once the orientation of the trimerization 
mode is set, the spontaneous polarization can emerge in only 
one direction, (d) Typical domain structure measured using 
piezoforce microscopy. The black and white regions corre- 
spond to opposite orientations of the ferroelectric polarization 
along the z axis. Note that the domain structure is isotropic, 
in spite of the layered crystal structure. 



tron microscopy to meet at vortex cores=s. Importantly 
(and unusually), while this mode lowers the symmetry to 
that of a polar space group it carries no net polarization, 
as any net local polarity vanishes macroscopically due 
to the non-zero mode wave vector. A secondary mode 
of symmetry (referring to the parent space group), 
which does not further lower the symmetry, provides the 
ferroelectric polarization (Fig. 1(c)). While the orien- 
tation of this secondary ferroelectric polarization is set 
by the in- or out- tilt of the K3 mode and does not re- 
sult in additional domains, it is essential for our exper- 
iments as it allows the straightforward imaging of the 
domain structure using piezoresponse force microscopy 
(PFM). Indeed PFM measurements reveal that domains 
of alternating polarization are locked to the trimeriza- 
tion domains around vortex core a^^'^^ yielding appealing 
six-fold patterns (Fig. 1(d)). Electric-field poling exper- 
iments have shown that the vortex cores are protected 
in the sense that they cannot be annihilated or driven 
out of the system by an electric field^ii^. Surprisingly, 
the domain structure and density of these topological de- 
fects when viewed from the side of the sample is similar 
to that viewed from the top in spite of the layered crys- 
tal structure and uniaxial ferroelectricity^ (Fig- 1(d)). 
This absence of anisotropy in the domain structure allows 
for straightforward determination of the defect densities 
from two-dimensional top-view scans of their areal den- 
sity, rather than requiring a complex three-dimensional 
analysis. 

First-principles calculations^ and Landau Theory 
analysis^S have shown that for small magnitudes of the 
trimerizing mode, the polar mode appears only as a 



FIG. 2: Mexican-hat potential energy surface of the hexagonal 
manganites. At high energy (the peak of the hat) the energy is 
independent of the angle of trimerization, and the system has 
(7(1) symmetry. At lower energy (in the brim of the hat), six 
of the trimerization angles become favorable (white circles), 
and the symmetry reduces to the sixfold discrete symmetry 
described by Zq. 



third-order term, and so the magnitude of the ferroelec- 
tric polarization just below Tc is vanishingly small. This 
is important for two reasons: First, the formation of the 
domain structure at Tc is not influenced by the system's 
attempts to minimize the depolarizing field from the fer- 
roelectric polarization. Strong evidence for this is given 
by the large numbers of electrostatically unfavorable 
head-to-head and tail-to-tail domain walls that form in 
i^MnOa , but rarely occur in conventional ferroelectrics^^. 
Second, first-principles calculations show that the energy 
lowering provided by the condensation of the K3 mode is 
independent of the angle of the tilt until the polar mode 
subsequently develops^22S. This means that the poten- 
tial below the phase transition temperature is given by 
the continuous "Mexican hat" form (Fig. 2). The dis- 
crete nature of the lattice does not manifest until lower 
temperatures when the domain structure is already de- 
termined. As a result we can use the mathematics of 
continuous symmetries which are usually assumed in the 
Kibble-Zurek mechanism. In this language, the full ro- 
tational symmetry is broken when the polyhedra tilt in 
the 2tt range of angles, resulting in a U{1) vacuum. 

For larger magnitudes of the K3 mode, obtained upon 
temperature decrease, a cross-over to linear coupling with 
the polar F^ mode occurs and the polarization becomes 
measurably large. This lifts the degeneracy with angle 
of the K3 mode and fixes the polyhedra into discrete 
tilt angles of 0,27r/3 or 47r/3, described by Z3 symme- 
try. The additional degeneracy provided by the direction 
( "in" or "out" ) of the polyhedral tilting gives an addi- 
tional Z2 symmetry reduction, resulting in Z2X = Zq. 
It is an open experimental question whether the onset 
of the F^ mode, which is observed ^ 300 K below Tc 
is an "emergence" or an additional isosymmetric phase 
transition2iiS2. 



3 



II. KIBBLE-ZUREK SCENARIO FOR RMNO3 

In this section we first show that the symmetry of 
i?Mn03 results in topologically protected vortex cores as 
described by the Kibble mechanism. We then analyze the 
vortex cores using the Zurek scenario to determine the 
density of topological defects that should be produced as 
a function of the cooling rate through the phase transi- 
tion. We use first-principles density-functional theory to 
evaluate the relevant parameters, and show that our pre- 
dictions are in excellent agreement with literature data. 



A. Kibble mechanism and the formation of 
topological defects. 

The requirements for the formation of topological de- 
fects at a phase transition within the Kibble mechanism^ 
are (i) a spontaneous symmetry breaking and (ii) a 
change in symmetry across the phase transition that cor- 
responds to a non-trivial homotopy group. The trimer- 
ization transition in i?Mn03 clearly fulfils the first con- 
dition; next we show that it also fulfils the second. 

As discussed earlier, in the temperature range just be- 
low the phase transition, iiMnOs exhibits a continuous 
symmetry; this allows us to use the methods and results 
of homotopy theory - which are developed for continu- 
ous symmetry groups - to assess the topology of i^MnOs. 
It is established within homotopy theory that the sym- 
metry characteristics of the order parameter, in our case 
U{1), can be used to assess the topological character- 
istics of a phase transition. To make the assessment, 
the order parameter symmetry is first mapped onto an 
ri-dimensional sphere. In the case of the U{1) symme- 
try, this is a one-dimensional circle, S^. Next a func- 
tion called the homotopy group, tt/j, which describes the 
topological nature of the order parameter symmetry is 
defined. If Wk differs from the identity, then it is non- 
trivial and topological defects are formed. It has been 
known since the 1960's^that nuiS^) is indeed non-trivial 
and in fact produces one-dimensional topological singu- 
larities, called strings or vortex cores^. Therefore the 
vortex cores in RMnOs are mathematically topologically 
protected, in concordance with their physical topological 
protection - their resistance to annihilation by an electric 
field - that we discussed earlie r^*^'^^ . We note also that 
within the Kibble mechanism the topological defects are 
remnants of the parent phase trapped within the lower 
symmetry phase. For i?Mn03 this implies that the high- 
symmetry paraelectric phase is preserved at the meeting 
point of the six domains defining a vortex core. 



B. Zurek scenario for iiMnOa. 

Within the Zurek scenario, the density of topological 
defects formed during a spontaneous symmetry breaking 



phase transition described by the Kibble mechanism fol- 
lows a power-law dependence on the rate at which the 
transition is crossed2i2^. In this section we first relate 
the material properties of iJMnOs to the parameters in 
the Zurek scenario. We then evaluate their magnitudes 
to predict quantitatively the temperature dependence of 
the defect formation within the Kibble-Zurek mechanism. 

Zurek's approach relies on the notion of competing 
timescales: The first relevant timescale is the time it 
takes for one region of the system to communicate its 
choice of vacuum state with another. This time becomes 
divergently long towards the critical temperature as the 
correlation length diverges, a phenomenon termed "crit- 
ical slowing down". The second is the quench time Tq 
that the system spends cooling through the phase tran- 
sition. The size of the domains is set at the tempera- 
ture Tc + ^Tf where the communication distance across 
which information can be transferred during the progress- 
ing phase transition becomes equal to the correlation 
length ^(T) (Fig. 3). While the correlation $(T) con- 
tinues to diverge as the temperature further approaches 
Tc, the communication length remains unchanged, and 
the system is unable to adapt to its increase. As conse- 
quence a "freeze-out" occurs in the temperature interval 
between Tc + AT/ and Tc - ATf. The size of the cor- 
related regions is unable to increase and so the domain 
size is fixed at the value ^{Tc + AT/) (Fig. 3). For fast 
cooling through the transition, the distance over which 
information can be transferred during the transition is 
small, and becomes equal to the correlation length at 
small values of ^(T). Therefore freeze-out occurs when 
the domain size is small (and consequently the number 
of topological defects is large). In contrast, for slow cool- 
ing, the distance for information transfer is large, and 
does not become equal to ^(T) until close to the phase 
transition temperature, where ^(T) is large. In this case 
large domains, with fewer topological defects, form. 

Here we summarize the derivation of the density of 
topological defects as a function of quench rate through 
the phase transition within the Kibble-Zurek mecha- 
nism. For a detailed derivation we recommend particu- 
larly Ref.— . First, we use critical scaling: As the system 
approaches Tc, the correlation length, ^, and relaxation 
time, r, diverge as 



r(T) = TO 



1-^ 

Tc 
T 

1 

Tc 



where ^0 is the zero-temperature correlation length and 
To is the zero-temperature time, which is equal to ^0 di- 
vided by the speed of information transfer in the system. 
Both ^0 and tq are system-dependent quantities, v and 
fi are critical exponents that are determined by the uni- 
versality class of the phase transition, that is its general 
behaviour as determined by the symmetry properties of 
the phase transition, irrespective of the material proper- 
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FIG. 3: Domain formation in the Kibble-Zurek scenario. 
Above Tc fluctuating regions of lateral extension d with uni- 
form orientation of the emerging order parameter occur (fuzzy 
patches). At high temperature (r> > Tc + AT/) the size of 
the correlated regions is determined by the correlation length 
(purple curve) . At temperature T = Tc + ATf a freeze-out of 
the lateral extension d begins. Below the freeze-out tempera- 
ture, Fluctuations continue to occur but the lateral extension 
of the fluctuating regions can no longer follow the diverging 
correlation length while the system continues to cool. The 
size of the fluctuation regions in this range, for example at 
r< remains set by the correlation length at the freeze-out 
temperature, (,{Tc + AT/), and corresponds to the "commu- 
nication length" which is the distance that information can 
propagate during the time in which the system cools from 
Tc + AT/ to Tc - M'f. Below Tc - ATf stable domains 
of the lateral extension ^(Tc + AT/) form (areas confined by 
black lines). 



ties of the specific system. 

Assuming that the temperature varies linearly with 
time near the phase transition, and taking t = at 
T — Tc, it is clear that T = Tc + t x r^, where Vq is 
the cooling rate. Rearranging this expression yields 
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T 



Tc T„ 



where t„ 



divided by the cooling rate, is called the 
"quench time" . 

The speed at which information is transferred in the 
material is then given by the characteristic velocity, 



c(T) 



t{T) 
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and the corresponding distance over which information 
can propagate in time t is 
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the last step we have substituted t 
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Equating the distance over which information can 
propagate to the correlation length yields an expression 
for the freeze-in temperature, Tf. 
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At temperature Tf, the domain sizes are frozen in with a 
characteristic length-scale given by the information prop- 
agation distance (which is equal to the correlation length) 
at the freeze-in temperature: 

0=Co(l + Ai-i^)^ i^^J 

Formally, this corresponds to the length-scale of the vor- 
tex "strings" , that is the three-dimensional continuation 
of the vortex cores through the bulk of the sample. The 
number of vortex intersections per unit area, n, is then 
approximately equal to the length of vortices per unit 
volume, 4-, giving 
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To apply the linear Zurek scaling law given in Eq. (jIIB[) 
to the hexagonal manganites we next identify the rel- 
evant time- and length-scales in the system, and eval- 
uate their magnitudes. Our electronic structure cal- 
culations were performed using density functional the- 
ory within the local density plus Hubbard- [/ approxi- 
mation following the Liechtenstein approach"^** with the 
double-counting corrections treated in the fully localized 
limit. Following previous literature studies^, we set the 
LDA+U parameters on the Mn 3d orbitals to U = 8 and 
J = 0.88 eV respectively and enforced an A-type antifer- 
romagnetic ordering. We used the projector-augmented 
wave method for core- valence partitioning"^^, which sig- 
nificantly reduces the required plane-wave energy cutoff, 
and carefully tested the convergence of plane- wave cutoff 
and k-point sampling. 

The zero-temperature correlation length, is usu- 
ally equated with the zero-temperature domain wall 
width in ferroelectrics. In order to extract this value 
we performed density functional calculations within the 
LBA+U method using the VASP code^'^'^'^. We con- 
structed supercells containing two 180° domain walls and 
in turn 120, 180, 240 and 300 atoms. We fixed the lattice 
constants of the supercells to those of the relaxed single- 
domain unit cell, and then optimized the internal posi- 
tions until the forces acting on all atoms converged to less 
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than 0.01 eV/A respectively. We initialized a different 
trimerization phase and ferroelectric orientation within 
adjacent domains, and subsequently perfomed full relax- 
ations on the structures; in all cases the system remained 
in the metastable multi-domain state. For all supercell 
sizes, we found that the structural phase defined by ei- 
ther the tilt of the MnO bipyramids or the direction of 
off-centring of the Y ions changes abruptly at the domain 
walls, indicating an effective domain wall width close to 
zero. This is consistent with a recent experimental elec- 
tron microscopy study—, and sets an upper limit on 
of - 1 A. 

To calculate the characteristic timescale of the system, 
To = — , we require s, which is the speed at which the 
system communicates the lattice distortion as it passes 
through the phase transition. For solid-state systems, s 
is given by the speed of sound. To calculate the speed 
of sound at zero kelvin we used the ABINIT^^ software 
package^SiiS to optimize the structure of a 10-atom unit 
cell, then calculated the full phonon band structure using 
frozen-phonon techniques. A supercell was constructed 
with doubling and trebling in each of the directions re- 
quired to sample the first Brillouin zone, then symmetry- 
distinct displacements were made to construct the full 
matrix of interatomic force constants. The dynamical 
matrix was diagonalized along each of the high symmetry 
lines shown in the phonon band structure using Fourier 
interpolation] "''^!"''^ The speed of sound was then extracted 
from the calculated phonon band structure by fitting the 
acoustic branch with a polynomial, then evaluating the 
group velocity 



dk 



k=Q 



which relates to 



through rq = 



— This is known 
with the exact number 



. The relevant velocity for our Kibble-Zurek fit is the 
doubly degenerate branch with the atoms displacing in 
plane and the wavevector propagating in plane. For this 
branch we obtained Vg = 640 m.s~^. 

Finally, for comparison with quenching experiments we 
also need the Curie temperature, Tq 
the cooling rate, 
experimentally to be ^ 1400 K 
depending on the R ion. 

Next we extract the critical exponents by identifying 
that, as in the case of the superfluid transition in *He, 
the i?Mn03 transition belongs to the universality class 
denoted as the 3D XY model. The values of the critical 
exponents for this universality class have been calculated 
using Monte Carlo simulations^"^ to be = 0.6717 and 
/i = 1.3132, giving a Kibble-Zurek scaling exponent of 
« 0.29. Taking the values introduced so far, with 
our upper limit for we find that domains of 5/im width 
should be formed for a quench time of Tq ^ 40 min (corre- 
sponding to a cooling rate of ~ 0.5 K/s), and domains of 
40 /im in around one month (cooling rate ~ 1.5 K/hour. 
These are readily experimentally accessible cooling rates. 

In Fig. 5(d) we compare our predicted scaling behav- 
ior with recently reported vortex densities measured as 



a function of cooling rate in ErMnOa (red circles)^. The 
agreement in scaling behavior between the experiment 
and the Kibble-Zurek prediction is excellent. Note in 
particular that the scaling exponent matches exactly the 
theoretical prediction; this is the first time to our knowl- 
edge that such good agreement has been obtained. We 
obtain the best match with a value of zero-temperature 
correlation length of 0.06 A (solid red line), consistent 
with the approximately zero domain wall width obtained 
in our density functional calculations. Hence, /JMnOa is 
the first system to show model Kibble-Zurek behavior in 
the laboratory. 

We thus find a unique situation in iiMnOs. Topologi- 
cally, it is a model system for the experimental verifica- 
tion of the Kibble-Zurek mechanism. In the temperature 
range where the Kibble-Zurek mechanism is expected to 
govern the formation of domains and the distribution of 
vortex-core singularities, the unwanted ferroelectric po- 
larization which could infiuence domain formation is ef- 
fectively absent. However, at room temperature the cou- 
pling of the distortive order parameter to the now finite 
ferroelectric polarization allows straightforward imaging 
of the topology and vortices via spatially resolved mea- 
surements of the ferroelectric domain structure. And fi- 
nally, a distinguishable range of domain sizes and hence 
defect densities is obtained for an experimentally accessi- 
ble range of cooling rates through the primary distortive 
phase transition at Tc- 



III. BEYOND THE KIBBLE-ZUREK LIMIT 

The Kibble-Zurek mechanism applies only to the 
regime in which the system has time to respond adia- 
batically to the cooling until the freeze-out temperature, 
Tc + ATf is reached. For faster quenching, it is expected 
that the Kibble-Zurek mechanism should break down and 
be replaced by a dynamics that is largely unknown. The 
unique combination of properties in i?Mn03 allows us to 
continue to increase the cooling rate beyond those ex- 
plored in Ref.'^ to investigate the evolution out of the 
Kibble-Zurek regime; we explore this range next. 

For our experiment we chose YMnOs rather than the 
ErMnOs that was used in Ref.— because of the greater 
thickness of our YMnOa samples. The experimental pro- 
cedure applied to the fiux-grown, c-oriented single crystal 
YMnOa platelets is shown in Fig. 4. For our anneal- 
ing experiments we used a conventional chamber furnace 
(carbolite CWF 13/13) which allows for temperatures 
from room temperature up to 1550 K exceeding consid- 
erably the Tc w 1270 K of our YMnOs samples. First 
we performed a prc-anncaling at 1420 K for 24 h under 
constant oxygen flow of 0.2 1/min. In the second step 
the samples were annealed again with a different cool- 
ing rate. Up to 8 K/min the temperature gradient was 
controlled by the furnace, higher rates were obtained by 
removing the fused-silica cell with the sample from the 
furnace and measured by an infrared camera (InfraTec 



6 



1 Pre-annealing 

1420 K 
24 h 

3 K/min/ \3 K/min 
O,,0.2 \lm\n\ 



2 Annealing at variable cooling rates 

1420 K 
24 h 

3 Klmin^ Variable cooling 

/O, , 0.2 l/minV\ rate x K/min 



^3 



3 Polishing 




4PFM 


Thin down by 
5 - 10 




Image domain 
distribution 





FIG. 4: Sketch of the anneahng procedure. After pre- 
anneahng the samples are heated to 1420 K with a hold 
time of 24 h and subsequent cooling at rates between 0.1 and 
1360 K/min. After the annealing cycle they are thinned down 
5-10 ^m by polishing in order to access the domain structure 
formed in the true three-dimensional bulk of the sample. 



VarioCAM) . The annealing was followed by thinning the 
samples down by 5-10 /zm using AI2O3 and polishing 
chemically-mechanically with 32 nm grain-size slurry. In 
the last step the ferroelectric domain structure was im- 
aged using piezoresonse force microscopy^'^®. We used 
a commercial scanning force microscope (Solaris, NT- 
MDT) operated in contact mode and applied an AC- 
voltage of 14Vpp at a frequency of ~40 kHz to a con- 
ductive Pt-Ir coated probe (NSC 35, Mikromasch) . The 
out-of-plane component of the piezoelectric response was 
recorded by the in-phase output channel of an external 
lock-in amplifier (SR830, Stanford Research) with a typ- 
ical sensitivity of 200 and time constant of 10 ms. 
Finally, the area density of vortices was extracted from 
the PFM images. 

In Fig. 5 (a-c) we show the measured distribution of 
domains in YMnOa for fast cooling rates of 3, 195, and 
1360 K/min, and in Fig. 5(d) we plot our measured 
vortex core densities as a function of cooling rate (blue 
circles). The independence of the vortex core density to 
the choice of R is confirmed by the smooth connection 
between the ErMnOs (red circles) and YMnOs (blue cir- 
cles) data points in Fig. 5. Following a crossover oc- 
curring between 0.3 and 3 K/min (see Fig. 5(d)) we 
observe a striking "anti-Kibble-Zurek" behavior in the 
measured density of domain vortex cores as a function of 
cooling rate: Within the fast-cooling regime, an increase 
in the cooling rate leads to a lowering of the density of 
vortex cores (larger domains); this behaviour is opposite 
to that predicted by the standard Kibble-Zurek mecha- 
nism. We emphasize that the behaviour is highly repro- 



ducible: We repeated our measurements on YMnOa sam- 
ples grown in different batches, and verified that there 
are no drift effects in consecutive annealing cycles. We 
therefore conclude that our results indicate an evolution 
out of the Kibble-Zurek regime at a cooling rate of ^ 
1 K/min in the hexagonal manganites. This crossover 
point corresponds to a correlation length (and hence a 
crossover domain size) of ^ 2.2 /im, and a relaxation 
time of ~ 4.1 x 10^"* s, with a characteristic informa- 
tion transfer velocity of ~ 5.4 x 10"'^ m/s, considerably 
reduced from the speed of sound by the critical slowing 
down. 



IV. DISCUSSION: POSSIBLE ORIGINS OF THE 
KIBBLE-ZUREK / "ANTI-KIBBLE-ZUREK" 
CROSSOVER 

A number of possible deviations from Kibble-Zurek be- 
havior have been discussed in the literature, but none 
of them are consistent with our measurements. Zurek^ 
showed that vortex-anti- vortex annihilation becomes sig- 
nificant at fast quench rates where domains are smaller 
and topological defects are closer together. Such vortex- 
anti-vortex annihiliation causes a leveling off of the rate 
at which the density of vortex cores increases with cooling 
rate, but not the decrease in density that we observe. The 
effect of nonlinearity in the quench rate on the density 
of defects has been calculated to yield a modified scaling 
law^, which again would not cause our observed turn- 
around at fast cooling rates. In Ref.- it was suggested 
that the observed production of defect-anti-defect pairs 
could be the result of a Kosterlitz-Thouless transition^, 
in which vortex-anti-vortex pairs are formed above the 
transition temperature and annihilate as the system is 
cooled. As a result, more vortex cores survive during a 
fast quench when the pairs do not have time to annihilate. 
This is the opposite of our observed fast-quenching be- 
havior. In addition, a Kosterlitz-Thouless system would 
show a dramatic change in the density of vortices after 
repeated annealing cycles as well as a dependence on the 
temperature at which the quench begins, neither of which 
we observe. 

A possible extrinsic infiuence on the domain struc- 
ture could be differences in chemical defect concentra- 
tion caused by the different cooling rates, such as off- 
stoichiometry, anti-site formation, or charge screening at 
the domain-walls. To test for this possibility, we heated 
our samples to within 2% of the the transition tempera- 
ture, and annealed them at this temperature for six hours 
under the conditions described above. No changes in the 
domain structure were observed on heating until 1270 K 
(just below Tc), at which point minor isolated domain 
wall movements, and, once, the formation of a vortex- 
anti-vortex pair were observed. On heating to 1320 K 
(just above Tc), a completely new, but statistically iden- 
tical domain pattern was obtained. These data and the 
aforementioned reproducibility of our data points in Fig. 
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FIG. 5: (a) to (c) Distribution of ferroelectric domains in z-oriented YMnOs samples after annealing cycles at different cooling 
rates. The images are obtained using PFM on an area of 30 x 30 fira^ and reveal a striking 'anti-Kibble-Zurek" behavior with 
higher cooling rates leading to larger domains, (d) Areal vortex core density as a function of cooling rate for slow cooling 
(red circles, Ref.— ) and fast cooling (blue circles, this work). Note the turnover in the cooling-rate dependence of the vortex- 
core density occurring between rates of 0.3 and 3 K/min. The red solid line is the result of our ab initio application of the 
Kibble-Zurek scenario with parameters from first-principles calculations. 



5 show that chemical drift effects do not play a role. 

We conclude therefore, that the most likely origin of 
the transition out of the Kibble-Zurek regime in i?Mn03 
is the breakdown of the assumption that the system is 
able to respond adiabatically to the cooling at high tem- 
peratures, above the freeze-out temperature Tc + ^Tf. 
Our observed behavior in the fast-cooling regime - slower 
cooling leading to a larger number of smaller domains - is 
of course reminiscent of nucleation-dominated behavior, 
with an activation energy for formation of the low sym- 
metry phase from the high symmetry phase. Nucleation- 
dominated phase transitions show charactistic first-order 
behavior, and a longer time spent at the transition al- 
lows a larger number of smaller domains to nucleate. We 
note that at the freeze-out temperature corresponding to 
the crossover quench-rate, the order parameter for the 

trimerization, rj ~ ^ "^tJ*^ ^ i has already reached 0.5% 
of its saturation value, taking the experimental value of 
ji = 0.29^. It is possible that this provides a sufficient 
discontinuity from the zero value to induce a first-order 
response. An alternative scenario is the fluctuation- 
induced first-order behaviour proposed for prototypical 
second-order phase transitions such as the normal-to- 
superconducting transition and for the nematic-smectic 
transition in liquid crystals*^, both of which belong to 
the same universality class - the 3D XY model - as the 
hexagonal manganites. Such an induced-first-order tran- 
sition could also explain the current controversy regard- 
ing the order of the trimerization transition in the hexag- 
onal manganites, with most experiments showing second- 
order behavior, but occasional reports of first-order char- 
acteristics. 



V. SUMMARY 

In Summary, we have shown that the multiferroic 
hexagonal manganites, -RMnOa, are model systems for 



testing the Kibble-Zurek scenario. Mathematically, they 
fulfil the symmetry requirements for the formation of 
topological defects, and practically, the defects are read- 
ily detectable, the quench rate can be varied over a wide 
range of relevant timescales, and extrinsic factors that 
might infiuence the phase transition behavior are absent. 
Our quantitative calculations of topological defect den- 
sity as a function of cooling rate using the conventional 
Kibble-Zurek model and parameters obtained using den- 
sity functional theory, yield excellent agreeement with 
literature data in the slow cooling limit where the conven- 
tional Kibble-Zurek mechanism is applicable. Our mea- 
surements of defect density at fast cooling rates, how- 
ever, reveal a surprising, apparently "anti-Kibble-Zurek" 
behavior in which faster cooling yields lower defect den- 
sities, reminiscent of a nucleation-dominated phase tran- 
sition. Since the expansion of the universe during the 
inflationary period was rapi d^^i^^ , the turn-around in 
vortex density which we flnd here at fast-cooling rates 
may be applicable, and the observed deficiency of cosmic 
strings in the cosmic microwave background'^ could be 
rationalised by the decrease in defect density observed in 
the anti-Kibble-Zurek regime. 
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